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We investigate, theoretically and numerically, the spin dynamics of a two-component exciton- 
polariton condensate created and sustained by non-resonant spin-polarized optical pumping of a 
semiconductor microcavity. Using the open-dissipative mean-field model, we show that the existence 
of well defined phase-locked steady states of the condensate may lead to efficient switching and 
control of spin (polarization) states with a non-resonant excitation. Spatially inhomogeneous pulsed 
excitations can cause symmetry breaking in the pseudo-spin structure of the condensate and lead 
to formation of non-trivial spin textures. Our model is universally applicable to two weakly coupled 
polariton condensates, and therefore can also describe the behaviour of condensate populations and 
phases in ’double-well’ type potentials. 
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I. INTRODUCTION 

Created in semiconductor microcavities, the two- 
component exciton-polariton condensates^ provide a 
rich playground for exploring the effects of internal 
(pseudo-spin) degrees of freedom on the dynamics of 
non-equilibrium superfluids. A remarkable range of spin- 
related phenomena has been explored in systems with a 
coherent resonant optical pump, which offers the possibil¬ 
ity to directly inject the spin state of exciton-polaritons 
and explore the ultrafast spin-switching dynamics^—. 
On the contrary, the non-resonant incoherent excitation 
responsible for the spontaneous formation of coherent 
exciton-polariton condensates is only beginning to be ex¬ 
perimentally explored in the context of spin-dependent 
effects. The reason for this disparity is a long-held belief 
that the incoherent reservoir (uncondensed hot exciton- 
polaritons), inevitably created by the incoherent exci¬ 
tations far above the energy of the condensed exciton- 
polariton quasiparticles, lacks spin-selectivity, i.e., owing 
to the fast spin-relaxation processes within the reservoir, 
the pump polarization cannot affect the spin of the con¬ 
densed polaritons. Thus, the theoretical modelling of an 
incoherently excited exciton-polariton system with spin 
degrees of freedom usually treats the reservoir as either 
fast£ or slowi response scalar density state, without any 
spin discrimination. 

Recent experimental studies of the optical spin-Hall 
effect^ employed microcavities with relatively short exci- 
ton lifetimes and succeeded in creating a spin-polarized 
reservoir of high-energy ’excitonic’ polaritons. Such a 


reservoir can be controllably replenished by a pump laser 
of the selected polarization, thus opening up the way 
to create and manipulate spinor condensates of exciton- 
polaritons in the regime of incoherent far off-resonant 
optical excitation^. 

In this work, we consider the exciton-polariton conden¬ 
sate in the regime of cw far off-resonant incoherent excita¬ 
tion, in the absence of external magnetic field. We anal¬ 
yse the non-stationary dynamics of condensate pseudo¬ 
spin and the relaxation dynamics towards possible steady 
states for both spatially homogeneous and inhomoge¬ 
neous cw pump. We show that the spin-polarized in¬ 
coherent reservoir allows efficient switching between dif¬ 
ferent spin states of the polariton condensate, both in the 
spatially homogeneous and some inhomogeneous pump¬ 
ing regimes. Furthermore, we analyse steady-state spin 
textures and demonstrate that excitation of steady states 
by a spatially inhomogeneous pulsed optical pump can 
lead to the creation of stable non-trivial spin textures. 

II. THE MODEL 

We model the spontaneously created polariton conden¬ 
sate by the open-dissipative Gross-Pitaevskii equations^ 
derived from the Hamiltonian for a quantum fluid with 
pseudo-spin degrees of freedom obtained, in different con¬ 
texts, in Refs^iii. The dynamical equations for the 
condensate wavefunctions corresponding to the left if- 
and right hand circular polarization states are writ¬ 
ten as follows (a = ±)* 
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Here m is the effective mass of the lower polarition, j c is 
the loss rate of polaritons which is the inverse of the po¬ 
larition life time, and R characterises the stimulated scat¬ 
tering rate from the reservoir into the condensate. The 
interaction constants u a and u b characterise the scatter¬ 
ing between the polaritons of the same and different po¬ 
larization (spin) states, respectively. It is well established 
that \u b \ <\u a ^. Also, g R characterises interactions be¬ 
tween the reservoir and condensate, where for simplicity, 
we assume that interactions between oppositely spin po¬ 
larized reservoir and condensed polaritons are negligible. 
The linear spin-mixing term proportional to J (Joseph- 
son coupling) accounts for the sample-specific linear po¬ 
larization splitting observed in experimental and usu¬ 
ally arising due to stress at the semiconductor interfaces. 

The equations for polariton wave functions are cou¬ 
pled to the equations for the reservoirs created by the 
incoherent cw pump with intensity P CT (r, t). It is usually 
assumed that all memory of the polarization of the inco¬ 
herent pump is losl— However, as mentioned above, 
it was shown that the creation of a spin-polarized exci- 
ton reservoir is possible with a spin-imprinting optical 
pump for sufficiently short-lived excitons. Because the 
precise spin-relaxation mechanisms for the reservoir are 
unknown, we will consider the simplest model of spin- 
dependent reservoir, which was introduced in— 

= P<r{r,t) — ('JR + R\4>a\ 2 ) n a , (2) 

where 7 # is the loss rate of the reservoir polaritons and 
P CT (r, t) is the pump. The dynamics of the reservoir is 
assumed to be fast on the scale of the polariton lifetime, 
i.e., comparable to the dynamics of the condensate, so 
that the model of slow-response, static reservoir— is no 
longer applicable. We note that the coherent resonant 
injection of spin— can be readily introduced into these 
equations. 

By introducing the characteristic time T = 7 — 1 , length 
L = ^h/(m"/ c ), and energy E u = Hjc scales^, we can 
rewrite Eq. CD and © in the dimensionless form: 

idti>cT = j-^V 2 +Ua\4’a\ 2 + UbIV’-o-l 2 + gma 

+ ^ l Rn <? - 7c] | 4>a + J^-a 
d t n a = P„(r) - ( 7 R + R\i/j a \ 2 )n cr 

where u a -P u a /(E u L 2 ), u b -> u h /{E u L 2 ), g R -A 
g R /{E u L 2 ), and J —> J/E u . In these scaling units, 
7 C = 1 , however, we will retain it in the subsequent for¬ 
mulas for clarity. In the following discussions, except 
some specified cases, all parameters take the numerical 
values indicated in Ref— -. 


III. HOMOGENEOUS STEADY STATES 

We start with the analysis of the spin-mixing dynamics 
for a spatially homogeneous condensate. In experiments, 
the quasi-homogeneous distribution of condensate den¬ 
sity can be achieved, e.g., by a flat-top ’super-Gaussian’ 
excitation. In this case, it is convenient to separate the 
amplitude and phase dynamics of the order parameter by 
using the transformation: ip a (t) = &a(t) exp ( i<j> a {t )) and 
introduce the density of the spinor components p a {t) = 
and the relative phase difference 9 = (f>- — </>+. Substi¬ 
tuting these spatially homogeneous states into Eq. ©^ 
we obtain the dynamical equations for the condensate 
density p a , the relative phase 9, and the reservoir den¬ 
sity n a : 

Pa = ( Rn a - "fc)Pa + cr2J^p a p-a sin 9 
9 = (u -{p + - p-) + g R 8n R (4) 

ha = Pa - (jR + Rpa)n a , 

where U = u a — u b and 6n R = n+ — n_. The dynamics 
described by Eq. © is rather universal, and would apply 
to the description of any two-mode polariton system with 
linear (Josephson-type) coupling such as condensates in a 
double-well potential^ or a two-level system^, provided 
that their spatial variation can be ignored. 

The amplitude and phase of each component in the 
condensate determine the polarization of the coherent 
photoluminescence emitted from the cavity. It is there¬ 
fore helpful to use the pseudo-spin representation (see, 
e.g.^L) which allows us to visualise the polarization of 
the condensate as a point on the Poincare (Bloch) sphere. 
The coordinates of the point are given by the components 
of the Stokes vector defined in the standard manner: 

S x — V P&P—& cos 9 , Sy — -\J paP—a sin 9 , 

1 1 ( 5 ) 

P+ - P -), •S'o = 2 (P+ + P-) 

with the total condensate density n c = p+ + p- defining 
the length of the Stokes vector: Sq = S 2 +S 2 +S 2 = n 2 / 4. 
Note that both the pseudo-spin and the total condensate 
density, S’o = n c / 2 , are functions of time, therefore it 
makes sense to consider the evolution of this vector in 
the normalized pseudo-spin space, s*(i) = Si(t) / So(t), 
i.e., on the Poincare sphere of unit radius. Although we 
will use the normalized Poincare sphere representation 
throughout the text, one should keep it in mind that its 
(not normalized) radius can change dynamically. 

The stationary states of the system m define fixed 
points on the Poincare sphere, which correspond to syn¬ 
chronized, phase-locked solutions of the Eq. © with a 
well-defined polarization state of the polaritons. Desyn¬ 
chronised states with a time-dependent phase correspond 
to Josephson oscillations, and, in the case of steady 
undamped oscillations to limit cycles on the Poincare 
sphere. 
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In principle, steady states can be obtained in a closed 
analytical form, however the expressions are cumber¬ 
some, and a more useful insight into the system’s dynam¬ 
ics is obtained in some limiting cases. As noted in£d£ ; 
the following conditions are of particular interest: (a) 
the excitons decay into the polariton channel dominates 
the relaxation dynamics of the reservoir, Rp± 7 _r; (b) 
the polariton interaction energy is greater than the tun¬ 
nelling energy, Un c J, which approximately defines 
the boundary of the synchronization region; (c) popu¬ 
lations of the two spin components are similar (semi- 
linearly polarized pumping), n c S z . Following the 
procedure outlined in£, one can derive the approximate 
self-consistent equation for the relative phase, which de¬ 
scribes a Josephson junction driven by a constant bias 
current (seei£, Ch. 8.5): 

+ y # = ( u + 9r)^y + UJn° c sme, (6) 

here 5P = P + — P-, n° c = (P + + P_)/ 7 C . Analy¬ 
sis of the phase space of this dynamical system shows 
that, for the equal pumping of the two spin compo¬ 
nents, i.e. linearly polarized pump, SP = 0, only phase- 
locked states exist. The system becomes desynchro¬ 
nised for an clliptically polarized pump when \I C \ = 
[\SP\/{P+ + P_)] [ 7 C (P + g R )/{2UJ)\ > 1, whereby the 
dynamics is dominated by limit cycles. For \I C \ < 1 the 
system admits fixed points determined by the condition 
sin 9 = I c and corresponding to the distinct steady states 
of the homogeneous system for each value of the pump 
intensity. 

Although Eq. © is obtained under a very restric¬ 
tive assumption of a nearly linear polarization of the 
polaritonic state, it can describe some experimental re¬ 
sults such as in Refj- 14 . Notably, when the pump power 
is large, conditions (a) and (b) can be automatically 
fulfilled once condition (c) is satisfied. Physically the 
condition (c) means that the system effectively evolves 
on a fixed-radius Poincare sphere, where its evolution 
is then given by the competition between the nonlin¬ 
ear polarition-polarition interaction and the Josephson 
tunnelling between two components. Meanwhile, if the 
pump power is large, then the polarition-polariton in¬ 
teraction dominates and Eq. © is applicable (detailed 
discussion can be found in Refji). 

In general, however, Eq. © does not always apply, es¬ 
pecially when the Josephson coupling dominates. There¬ 
fore, in what follows we will consider steady states of Eq. 
© given by the pumping conditions of different polarza- 
tion and strength. 

A. Linearly polarized pump 

Most commonly, the optical pump producing a polari¬ 
ton condensate in non-spin-resolved experiments is lin¬ 
early polarized. In our model, this means that the pump 
is balanced: P + = P_. Under this condition, from the 



FIG. 1: (Color online), (a)-(c) Fixed point solutions of Eq. (J4|) 
under balanced (linearly polarized) pumping. Branch 1 and 3 
are elliptically polarized solutions; branch 2 and 4 are bond¬ 
ing (si) and anti-bonding (S 2 ) solutions respectively. Only 
the solid branch is sable (see text). Green dots: correspond¬ 
ing points to the inhomogeneous simulation, (d) (see Sec. IIVI) 
Inhomogeneous simulation results of Eq. ©• Evolution of the 
integrated Stokes parameters in phase space under balanced 
Gaussian-shape pumps. Parameters (see Sec. IIVI for the defi¬ 
nitions): P = 5, = 0.5, a x = a y = 5, and J=0.5. Red dot: 

initial state of evolution; green dot: final state of evolution. 
Inset: energy evolution. 

first of Eq. ©, it follows that the synchronized state 
can be reached when Rn a = 7 C and the internal Joseph¬ 
son current Ij = 2J yjp+p- sinf? becomes zero. Conse¬ 
quently, the phases are locked to 9 = 0 or 9 = n (mod 
2n). The two fixed points on the Poincare sphere corre¬ 
sponding to these phase locked linearly polarized states 
states are Si = (1,0,0) (corresponding to 9 = 0), and 
s 2 = (—1,0,0) (corresponding to 9 — 7r). The density 
of the steady state grows linearly with the pump rate: 
Pa = {Pth/ 7c) {Pa/Pth ~ 1), see Fig. □ (a-c), where the 
normalized pump power is defined as P a = P a /Pth and 
for the linearly polarized pump particularly P + = P_ = 
P. 

Following the standard Lyapunov stability analysis^, 
we deduce that, for our choice of the sign of J, the n 
out-of-phase state S 2 (the so called anti-bonding state) 
is stable and the in-phase state si is unstable. Direct 
numerical integration of Eq. © with different initial 
conditions also confirm that S 2 is stable Fig. ©a-b). We 
note that this selection of the anti-bonding steady state 
in a weakly (linearly) coupled two-state system has also 
been recently confirmed i n 18 i 19 . This effect persists even 
for the inhomogeneous (Gaussian) pumpingiS, as shown 
in Fig. Md) and discussed in Sec. m below. 
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FIG. 2: (Color online). Trajectories of the pseudo-spin vec¬ 
tor obtained from numerical solutions of Eq. Q with J = 0.5 
and various interaction strengths between the condensate and 
reservoir gn. Red dots: initial states of evolution; Green 
dots: final states of evolution. Initial conditions are arbitrary. 
Parameters are: (a )gn = 1.5 x 1O _2 ,0(O) = 0.5; (b)t/.R = 
1.5 x 10 -2 ,0(0) = 2; (c) Artificially enlarged gn = 1.5 x 10 _1 ; 
(d) Artificially reduced gu = 1.5 x 10 -3 . 


The other steady states arise when p a > 
{Pth he) {Pa/Pth - 1), where the sign of the inequality 
is the opposite for the other polarization component, see 
Fig. □ (a-c) dashed branch 1 and 3. They are sustained 
by a non-zero internal Josephson current which will grow 
with the pumping power. Lyapunov analysis reveals 
that they are unstable. 

Introduced by the interaction between the condensate 
and the reservoir, the blue shift term which is propor¬ 
tional to gu plays an important role in the spin dynamics. 
It can modify the fixed point properties when its energy is 
compatible with the Josephson tunnelling energy. Start¬ 
ing with some arbitrarily chosen small amplitude initial 
conditions, Fig. EJc-d) illustrates how different values of 
gu can affect the evolution of Stokes parameters in phase 
space. Indeed, both the positions of fixed points and 
their stability can change, e.g., when g r = 1.5 x 10 , 

both s 2 and si become stable, see Fig. Eld); higher value 
of gn will lead to the system in an elliptically polarized 
state, see Fig. He). Here, the trajectories are obtained 
by numerical integration of Eqs. ©• 


B. Elliptically polarized pump 

When the imbalanced polarization-selective pumping 
is introduced, the dependence of the density of the or¬ 
thogonal spin components on the pump power displays a 
hysteresis-like profile. However, it does not correspond to 
the typical bistability behaviour, since most of the (light 
blue) curve corresponds to an oscillatory unstable fixed 
point. Fig. [3] (a-c) shows a typical fixed point distribu- 






FIG. 3: (Color online), (a)-(c) Selected fixed point solutions 
of Eq. ((4]), with an elliptically polarized pump for P+ = 2.4. 
Curves with darker color are stable while those with lighter 
color are unstable. Red (green) dot: the steady state be¬ 
fore (after) a pulse is added, (d) Inhomogeneous simula¬ 
tion results of Eq. © (see Sec. m Evolution of integrated 
Stokes parameters in phase space. Pumping parameters: 
P = 4, Ph = 0.4, a x = a y = 5, J = 0.5; Pulse parame¬ 
ters: P se = 0.12, Ph = 1 ,a px = a py = 5. (see Sec. IIVI for 
the definitions.) Red dot: initial state of evolution; green dot: 
final state of evolution. Inset: energy evolution. 


tion against small pumping imbalance, where all unsta¬ 
ble solution branches with at least one vanishing p a are 
not shown for clarity. Stable fixed points exist in semi¬ 
circular polarized pumping region (dark blue), which will 
be discussed in the next subsection, and the semi-linearly 
polarized pumping region (dark red), where at the point 
P_ = P + the fixed point coincides with the single anti¬ 
bonding stable state p + = p- shown in Fig. [TJ 

One should note that the above discussion is based 
on the Lyapunov stability analysis, which describes the 
linear stability to a long-wavelength (spatially homoge¬ 
neous) perturbation with the wave vector k = 0. The 
steady state may become unstable to a spatially mod¬ 
ulated perturbation with k ^ 0, see, e.g., discussions 
about the modulational instability in 2D£ and 1D22 cases. 

As a verification of the existence of the stable branches, 
we performed a pulsed excitation of a steady state in Fig. 
[3] (a-c) where the red and green dots represent stable 
steady states before and after the pulse. They fit the in¬ 
homogeneous simulation result, Fig. [3] (d), qualitatively 
well, (see Sec. |TV]for details.) 




















IV. INHOMOGENEOUS SPIN-SWITCHING 
DYNAMICS 
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FIG. 4: (Color online). (a)-(c) Fixed point solutions of 
Eq. (|4)l. with P_ = 0. Only the solid branch is stable (see 
text). Green dots: corresponding points to the inhomoge¬ 
neous simulation, (d) (see Sec. lIVI) Inhomogeneous simulation 
results of Eq. ©. Evolution of integrated Stokes parameters 
in phase space under circularly polarized pump. Parameters: 
P = 3, Pi m = 0, a x — a y = 5, and J=0.5. Red dot: initial 
state of evolution; green dot: final state of evolution. Inset: 
energy evolution, (see Sec. IIVI for the definitions.) 


C. Circularly polarized pump 


The steady state shown in Fig. [3] (a-c) at P_ = 0 
corresponds to the situation when only one spin compo¬ 
nent is pumped, which is physically akin to the strongly 
imbalanced dissipative double-well Josephson junction!!. 
(The opposite case where T+ = 0 is similar.) In this lim¬ 
iting case, the fixed point solutions of Eq. are shown 
in Fig.[4l where the steady state for P + = 2.4 is the same 
as that represented by the dark blue curve in Fig. [3] (a-c) 
at P- = 0. 

The stable branch whose 9 is small but never van¬ 
ishes corresponds to the so called self-trapped state with 
a strongly imbalanced population. It is a consequence 
of nonlinear interaction s 21 : 22 and, in the context of a 
double-well polariton system, has been studied both 
experimentally!! and theoretically!!. 

Again, Fig. [I] (d) shows us the inhomogeneous simula¬ 
tion result verifying the position of fixed points predicted 
by the homogeneous model. 


In the spatially inhomogeneous situation, it is well 
known that the in-plane magnetic fields, either being ex¬ 
ternally applied or caused by the TM-TE splitting!!, can 
lead to nontrivial spatial structures for both density and 
spin distributions!!!!!. The linear spin coupling can also 
be regarded as an effective magnetic field pointing along 
the x axis, which can lead to the appearance of spin pat¬ 
terns. However, here we analyse steady states of the 
polariton spinor systems sustained by a cw pump with 
a Gaussian shape and show that the spatially averaged 
evolution of the spin components qualitatively follows 
the prediction of the homogeneous model discussed in 
Sec. mi The results obtained within the homogeneous 
approximation can therefore be used as a guide for spin¬ 
switching manipulation with a realistic, spatially inho¬ 
mogeneous pump. 

To determine the evolution of pseudo-spin under spa¬ 
tially inhomogeneous excitation, we solve directly the 
open-dissipative GP equations d3j) with a radially sym¬ 
metric Gaussian pump by using the split-step method!!. 
The total intensity of the pump is defined as P(x,y) = 
P 0 e (x ' a x +y ! a y>, where a x<y is the width along the x 
and the y direction. The pump for each circularly po¬ 
larized component is given by splitting the energy from 
P(x,y), which is P(x,y) = P + (x,y) + P-(x,y), since in 
many experiments the pump is given by a single laser—. It 
is convenient to use the spin bias defined as P,/ m = Pff / Pq 
to denote those two pumps, then P_(x,y) = P// n P(x,y) 
and P + (x,y ) = (1 — P-f n )P(x,y). Also, throughout the 
text, the notation of the normalized (inhomogeneous) 
pump is given by, P = Pq/P t h (normalization against its 
peak value), where the pump threshold P t h is the same 
as that in Sec lIIIl 

By adjusting the pumping power and the imbalance 
value, with an arbitrary small amplitude initial condi¬ 
tion, the system can be driven to a synchronized steady 
state without any external potential. The condensate 
mean-field energy functional, E = E + + £L, serving 
as the measure of kinetic and interaction energy of the 
condensate!!, can be used to detect whether the steady 
state is reached: 

Ea = j d r { \V^ + u M * + 9 R n„\U + ^ 

+ u b \ijj a \ 2 \ip - a \ 2 + JRe(ip*ip a )\, 

where the integration is over the area where the conden¬ 
sate density is non-negligible. 

To map the evolution of the inhomogeneous conden¬ 
sate to that of the Stokes parameters, we trace the evo¬ 
lution of the spatially integrated quantities as s - nt (f) = 
/ Si(x,y)dr/S q 1 * , where Sq 1 * is given by the spatially 
integrated length of the Stokes vector. The integration 
area is the same as that used to determine the steady 
state energy ©■ In the following, unless it is indicated 
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FIG. 5: (Color online), (a)-(c) Synchronized steady state, 
with P = 4, P- zn = 0.4, a x = a y = 5, and J = 0.5: (a) Density 
profiles (along y = 0) and spatial distribution of s z for the 
steady state, (b) Evolution of integrated Stokes parameters 
and energy, (c) Evolution of integrated Stokes parameters in 
phase space. Starting form the red dot; ending with the green 
dot. (d) desynchronised state, with P = 8, P^ = 0.4, a x = 
a y = 5, and J = 0.5: evolution of integrated Stokes param¬ 
eters. Starting form the red dot; ending with a closed or¬ 
bit.Inset: energy evolution. The integrations were over areas 
where n c (x , y) > 10~ 3 . 


explicitly; we shall refer to the normalized s* nt (f) sim¬ 
ply as integrated Stokes parameters. This normalization 
procedure, although eliminating information on the spa¬ 
tial distribution of the pseudo-spin, allows us to recover 
the spatially averaged polarization state. This averaged 
spin dynamics could be observable experimentally even 
without performing spatially resolved polarinretry of the 
cavity photoluminescence. 

Polarization dynamics for the linearly and circularly 
polarized pump are shown in Fig. [T^d) and Fig. [2(d), 
respectively. As we can see, the spatially averaged dy¬ 
namics in the case of inhomogeneous pumping fits the 
homogeneous predictions well. 

For a typical steady state under elliptically polarized 
pumping, the value of the cross section densities and spa¬ 
tial distribution of the pseudo-spin is shown in Fig. [5](a). 
Its build up process is similar to that of the linearly- 
polarized pumping case, see Fig. ©T-c). The linear cou¬ 
pling J causes switching of the dominant density between 
the two components as polaritons spread out from the 
pump source, which leads to the appearance of radially 
symmetric domains of polarization density s z . The larger 
J leads to the denser spin pattern and loss of correspon¬ 
dence between the averaged spin dynamics in inhomoge¬ 
neous and homogeneous cases. If J is sufficiently small, 
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FIG. 6: (Color online). Steady state integrated quantities 
against the pump imbalance P F,. Parameters: P = 4, a x = 
a y = 5, J = 0.5. Between the gap 0.42 < P^ m < 0.58, one 
needs to increase the P value to reach steady states (see text). 


the polarization state is almost homogeneous across most 
of the integration area. 

The degree of control of the averaged polarization state 
of the condensate attainable by an incoherent, far off- 
resonant spin-polarized excitation according to our model 
is presented in Fig. [G] and could be tested in an exper¬ 
iment, which would validate our assumptions. Due to 
our choice of the pump power distribution configuration, 
in Fig. © between 0.42 < P^ rl < 0.58 one needs to in¬ 
crease P in order to reach steady states. With a larger 
pump power, however, the degree of control of polariza¬ 
tion decreases because the system will easily fall into a 
desynchronised state even with small value of imbalance. 

For example, for an elliptically polarised pump, within 
the region marked with ’D’ (which stands for desynchro¬ 
nised) in Fig. © none of the fixed point branches are sta¬ 
ble, which does not rule out the existence of the closed 
orbits (limited cycles) on the Poincare sphere correspond¬ 
ing to desynchronised states. As shown in RefJ^, for a 
desynchronised state the density p a (t) maintains peri¬ 
odic oscillations and the evolution of its Stokes vector 
is similar to Fig. [5] (d) where the angle of the trajec¬ 
tory plane can vary depending on the system parame¬ 
ters. This regime cannot be accessed by Eq. © because 
the latter can only describe limit cycles in the vicinity of 
the s z = 0 plane. Indeed, it can be revealed by inhomo¬ 
geneous simulations of the full model equations that, in 
such a state, the spatial distribution of the spin pattern 
keeps breathing. The corresponding Stokes parameters 
circle around the Poincare sphere^, reducing the degree 
of (time-averaged) polarization of its luminescence. A 
similar pinning and depinning effect has been observed 
experimentally^. Here this effect is associated with mov¬ 
ing in and out of the phase-locked synchronized regime of 
the condensate dynamics caused by intrinsic interactions 
between the two spin components, so that it can arise 
even without taking into account the structural disorder 
in the sample. 

The coexistence of stable and unstable fixed points 
described in Sec. m indicates that an inhomogeneous 
pump might be used to switch the system between 
different density and polarization states. Indeed, we 
find that the efficient switching between different spin 
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FIG. 7: (Color online). Spin switching dynamics for the case 
of Gaussian pump, (a-c) Energy of the initial and final state 
after the switching pulse. Dash line: pulse switched on. Dash- 
dot line: pulse switched off; (d-f) Integrated total density and 
integrated Stokes parameters during the switching process for 
switching between the steady states shown in (c). Parameters: 
the same as Fig. [3] (d), except for (a) P se = 0.4, P]^ = 0.65; 
and (b) P se = 0.9, P.' m = 0.65. 


and density states can be achieved by applying a radi¬ 
ally symmetric incoherent pulsed excitation of the form: 
-P S e(r, t) = (l/4)P 0 se {l+tanh[r p (t-to)]}{l + tanh[r p (ti- 
f)]}e - ( x / a p* +v / a py\ where Pq c is the peak value, a px , py 
determines the pulse width, to and t\ (t\ > to) are the 
pulse switch on and switch off time respectively, and t p 
is a coefficient controlling adiabaticity of the excitation 
(large t p corresponds to a non-adiabatic pulse). The nor¬ 
malized notation for the pulse is given by P se = Pq £ /P th- 
Note that the combined intensity after the pulse was 
added is given by the interference between pump and 
pulse modes. Unless specified, they are regarded as be¬ 
ing in phase. A radially symmetric non-adiabatic pulse, 
a px = CL py , causes an abrupt change in the condensate en¬ 
ergy and a strong outward density flows combined with 
the internal Josephson currents leading to spin mixing 
dynamics [Fig. CO) d-f)]. As a result, the system can enter 
a different steady state [Fig. [7](a,c)], or a phase desyn¬ 
chronised state [Fig. [7(b)] . If the original steady state 
is stable, the perturbed system restores its initial density 
and polarization state after the pulse is switched off. 

Regardless of the fact that with strong spatial vari¬ 
ation the homogeneous results would not be generally 
applicable, one can still compare the spatially-averaged 
inhomogeneous results of Fig. [7(c-f ) with the fixed point 
solutions of Eq. as long as the former begins and 
ends with steady states. Fig. [3] shows two of the rele¬ 
vant fixed point solution branches that are qualitatively 
comparable to the simulation data. The initial (red) and 
final (green) steady states correspond to the flipping of 
s x shown in Fig. Etd) before and after the pulse was 
added. These results are ready to be tested in experi- 



FIG. 8: (Color online). Spin wave excitation for the case 
of an elongated pump with the imbalance Pt rn = 0.4 and 
a tightly focused radially symmetric pulse with P.jp„ = 0.65 
and different widths (see text), (a) Density profile of the ip+ 
component showing a density dip associated with a phase fold 
(b). (c) Cross-section of density and phase across the fold, (d) 
Density profile of the ip- component corresponding to (a); 
(e) Density profile of the ip+ component showing stationary 
configuration of vortex pairs connected by a phase fold (f); 

(g) Cross-section of density and phase across the phase fold. 

(h) Density profile of the ip- component corresponding to (e). 
All density plots are logarithmic. 


ments, while the value of pump power might be given by 
a spatially averaged one depending on the specific shape 
of the pumping laser. 


V. EXCITATION OF NON-TRIVIAL SPIN 
TEXTURES 


While the density and spin wave excitation in the 
regime of a radially-symmetric cw pump and pulse leads 
to the efficient switching between different spin and den¬ 
sity states, the spatial dynamics can be captured by the 
homogeneous or spatial averaging approximation, as long 

















































as the radial symmetry of the spatial density distribu¬ 
tions for both spin components is preserved. However, 
the situation changes dramatically when the radial sym¬ 
metry of the system is broken. This can be achieved by 
combining a radially symmetric normal incidence pump 
with elongated excitation, e.g., at a steep angle to the 
sample surface, or vice versa. Fig. [8] demonstrates the 
consequence of pulsed perturbation of a steady state es¬ 
tablished with an elongated (a x = 16, a y = 5) pump. 
The pulse is a strong, off-resonant, tightly focused Gaus¬ 
sian beam with a px = a py = 4 (a-d) or a px = a py = 3 
(e-h) and peak power comparable to that of the cw pump. 
It was set to having 7r/2 phase difference with the pump 
and thus no interference pattern appeared. Remarkably, 
the symmetry of the density flows caused by the pulse 
is now broken, leading to different flow speeds along the 
two symmetry axes. As a result, larger area pulse pro¬ 
duces a phase fold associated with a pronounced dip in 
the condensate density Fig. [5] (a-c). An even tighter 
pulse will result in the two phase folds terminating on 
vortices, leading to formation of a stable and stationary 
configuration of two vortex-antivortex pairs. 

Although nontrivial phase structures appear only in 
one of the spin components [see Figs. |8](d), (h)], they are 
maintained and stabilised by internal Josephson current 
between the spin components and do not exist in the 
absence of linear coupling (J = 0). We note that no phase 
singularities appear in the linearly-polarized basis, which 
sets these structures apart from the previously described 
half-charge vortices22r— 

VI. CONCLUSIONS 

In conclusion, we have investigated in detail the forma¬ 
tion and polarization structure of the steady states of an 
exciton-polariton condensate created by far off-resonant 
spin-polarized pump in a semiconductor microcavity with 
linear polarization splitting. We model the dynamics 


of the system under the assumption that the incoher¬ 
ent pump creates a fast-responding and rapidly decay¬ 
ing reservoir, which significantly affects the condensate 
dynamics. Our analysis includes effectively spatially ho¬ 
mogeneous pumping conditions, as well as spatially in¬ 
homogeneous pumping with a Gaussian profile. The po¬ 
larization dynamics in the latter case can be mapped to 
the homogeneous dynamics through a spatially averag¬ 
ing procedure. We have demonstrated that the phase¬ 
locking conditions for existence of stationary states nat¬ 
urally leads to the self-trapped states with a nontrivial 
phase relationship between the condensate components, 
including 7r out-of-phase (’anti-bonding’) states. The ex¬ 
istence of both phase-locked and spin-beating states can 
enable efficient control and switching of the polarization 
states in the system. Experimentally, these findings can 
be tested with spatially averaged polarization measure¬ 
ments under incoherent spin-selective excitation condi¬ 
tions. 

In addition, we have analysed formation of non-trivial 
spin-textures under pulsed, spatially inhomogeneous ex¬ 
citation, and demonstrated numerically reliable creation 
of stable vortex-antivortex pairs. 

Our analysis, with suitable modification of parameters, 
is widely applicable to a general two-state polariton sys¬ 
tem with weak linear coupling, including weakly linked 
spatially separated condensates or multi-mode conden¬ 
sates in shallow potential traps. 
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